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Abstract 

A first measurement is presented of exclusive photoproduction of p° mesons associated 
with leading neutrons at HERA. The data were taken with the HI detector in the years 
2006 and 2007 at a centre-of-mass energy of y/s = 319 GeV and correspond to an inte¬ 
grated luminosity of 1.16 pb 1 . The p° mesons with transverse momenta p-j < 1 GeV 
are reconstructed from their decays to charged pions, while leading neutrons carrying a 
large fraction of the incoming proton momentum, xl > 0.35, are detected in the Forward 
Neutron Calorimeter. The phase space of the measurement is defined by the photon vir¬ 
tuality Q 2 < 2 GeV 2 , the total energy of the photon-proton system 20 < W 1P < 100 
GeV and the polar angle of the leading neutron @n ^ 0 .75 mrad. The cross section of 
the reaction yp —» p°mr + is measured as a function of several variables. The data are 
interpreted in terms of a double peripheral process, involving pion exchange at the pro¬ 
ton vertex followed by elastic photoproduction of a p° meson on the virtual pion. In the 
framework of one-pion-exchange dominance the elastic cross section of photon-pion scat¬ 
tering, cr el (77r + —» p u n + ), is extracted. The value of this cross section indicates significant 
absorptive corrections for the exclusive reaction yp —>• p°mr + . 
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1 Introduction 


Measurements of leading baryon production in high energy particle collisions, i.e. the produc¬ 
tion of protons and neutrons at very small polar angles with respect to the initial hadron beam 
direction (forward direction), are important inputs for the theoretical understanding of strong 
interactions in the soft, non-perturbative regime. In ep collisions at HERA, a hard scale may be 
present in such reactions if the photon virtuality, Q 2 , is large, or if objects with high transverse 
momenta, p T , are produced in addition to the leading baryon. In such cases the process usually 
can be factorised into short-distance and long-distance phenomena and perturbative QCD often 
is applicable for the description of the hard part of the process. 

Previous HERA measurements [QH71 have demonstrated that in the semi-inclusive reaction 
e + p —> e + n + X the production of neutrons carrying a large fraction of the proton beam 
energy is dominated by the pion exchange process. In this picture a virtual photon, emitted 
from the beam electron, interacts with a pion from the proton cloud, thus giving access to the 
7*7r cross section and, in the deep-inelastic scattering regime, to the pion structure function. 

The aim of the present analysis is to measure exclusive p° production on virtual pions in 
the photoproduction regime at HERA and to extract the quasi-elastic 77r —>■ p°7r cross section 
for the first time. Since no hard scale is present, a phenomenological approach, such as Regge 
theory |JS), is most appropriate to describe the reaction. In the Regge framework such events 
are explained by the diagram shown in figure [TJi which involves an exchange of two Regge 
trajectories in the process 2 —> 3, known as a Double Peripheral Process (DPP), or Double- 
Regge-pole exchange reaction This process can also be seen as a proton dissociating into 
(n, 7r + ) system which scatters elastically on the p° via the exchange of the Regge trajectory with 
the vacuum quantum numbers, called the “Pomeron”. 

In the past, similar reactions were studied at lower energies in nucleon-nucleon and meson- 
nucleon collisions [HO]4l4l . Most of the experimental properties of these reactions were suc¬ 
cessfully explained by the generalised Drell-Hiida-Deck model (DHD) [[1514171 . in which in 
addition to the pion exchange (figure [TJi) two further contributions (figure [Tb. [[}.') are included. 
The graphs depicted in figures Qi> and [lk give contributions to the total scattering amplitude 
with similar magnitude but opposite sign lfT8llT9l . Therefore they largely cancel in most of the 
phase space, in particular at small momentum transfer squared at the proton vertex, t —> 0, such 
that the pion exchange diagram dominates the cross section |[T71 . One of the specific features 
observed in these experiments is a characteristic t' dependence at the ‘elastic’ vcrtcxQ. with 
the slope dependent on the mass of the (nir) system produced at the other, pnir + , vertex, and 
changing in a wide range of approximately 4 < b(m) < 22 GcV 2 . The Deck model in its 
original formulation cannot fully describe such a strong mass-slope correlation and interference 
between the amplitudes corresponding to the first three graphs in figure Q] has to be taken into 
account to explain the experimental data ||20[1211 . 

In the analysis presented here only the two charged pions from the p° decay and the lead¬ 
ing neutron are observed directly. The pion from the proton vertex is emitted under very small 
angles with respect to the proton beam and escapes detection. This leads to a background con¬ 
tamination from events with a different final state, which originate from diffractive dissociation 

'in the present analysis elastic vertex corresponds to the p°!P vertex, figureQ] 
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Figure 1: Generic diagrams for processes contributing to exclusive photoproduction of p° 
mesons associated with leading neutrons at HERA. The signal corresponds to the Drell-Hiida- 
Deck model graphs for the pion exchange (a), neutron exchange (b) and direct pole (c). Diffrac¬ 
tive scattering in which a neutron may be produced as a part of the proton dissociation system, 
M Y , contributes as background (d). The N* in (c) denotes both resonant (via A r+ ) and possible 
non-resonant n + n + production. 


of the proton into a system Y containing a neutron (figure [Hi). Using the HI detector capabil¬ 
ities in the forward region such processes can be suppressed to a certain extent. The residual 
background contribution is estimated from a Monte Carlo model tuned to describe vector meson 
production in diffractive dissociation at HERA. 

The analysis is based on a data sample corresponding to an integrated luminosity of 1.16 
pb 1 collected with the HI detector in the years 2006 and 2007. During this period HERA col¬ 
lided positrons and protons with energies of E e = 27.6 GeV and E p = 920 GeV, respectively, 
corresponding to a centre-of-mass energy of y/s = 319 GeV. The photon virtuality is limited 
to Q 2 < 2 GeV 2 with an average value of 0.04 GeV 2 . 


2 Cross Sections Definitions 

The kinematics of the process 

e(k) + p(P) -A e(k') + p°{V) + n(N) + vr+, (1) 

where the symbols in parentheses denote the four-momenta of the corresponding particles, is 
described by the following invariants: 

• the square of the ep centre-of-mass energy s = (P + k) 2 , 

• the modulus of the four-momentum transfer squared at the lepton vertex 

Q 2 = -q 2 = ~{k - k r ) 2 , 

• the inelasticity y = (q ■ P)/(k ■ P), 

• the square of the yp centre-of-mass energy W 2 p = (q + P) 2 ~ ys — Q 2 , 

• the fraction of the incoming proton beam energy carried by the leading neutron 

x L = (q- N)/(q ■ P) ~ E n /E p , 
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• the four-momentum transfer squared at the proton vertex 

f = (p _ ft )2 ^ ^ 

' ' XL XL ’ 

• the four-momentum transfer squared at the photon vertex t' = (q — V) 2 . 


Here E p: m p ,E n , m n represent the energy and the mass of the incoming proton and the outgoing 
leading neutron, respectively. The square of the y7T centre-of-mass energy is then given by 

WP^W* p (l-x L ). 

Experimentally, the kinematic variables at the photon vertex (the mass M p , the pseudora¬ 
pidity ? 7 P and the transverse momentum squared p\ p of the p° meson) are determined from the 
p° decay pions, while those at the proton vertex (x L and Pr TI ) are deduced from the measured 
energy and scattering angle of the leading neutron. 

In the limit of photoproduction, i.e. Q 2 —> 0, the beam positron is scattered at small angles 
and escapes detection. In this regime the square of the pp centre-of-mass energy can be recon¬ 
structed via the variable W 2 = s y rec , where y rec is the reconstructed inelasticity, measured 
as y rec = (. E p — p Z)P ) / (2E e ). Here, E p and p z>p denote the reconstructed energy and the mo¬ 
mentum along the proton beam direction (z-axis) of the p° meson and E e is the positron beam 
energy. The variable t' can be estimated from the transverse momentum of the p° meson in the 
laboratory frame via the observable t' rec = —p 2 rp to a very good approximation^. 

The cross section of the exclusive reaction © can be expressed as a product of a virtual 
photon flux / 7 / e and a photon-proton cross section a rp : 

dydQ 2 = ) (T w( W ip(y))- (2) 

In the Vector Dominance model (VDM) 112211231 taking into account both transversely and lon¬ 
gitudinally polarised virtual photons the effective photon flux is given by 


fy/eiV, Q 2 ) 


a 

2irQ 2 y 


i + (i -yf 


2(1 - y) 




(3) 


where a is the fine structure constant and in = m 2 e y 2 /(1 — y), with m e being the mass of the 
electron and M p is the p° meson mass. 

In the one-pion-exchange (OPE) approximation ll24ll . which is valid for very small pj n ^ 
m 2 , the photon-proton cross section can be further decomposed into a pion flux, describing 
p —y nir + splitting, convoluted with a photon-pion cross section: 


d (lldyp, XX., t ) 
dxLdt 


fn/p{,-^Li t') (T 77 r (Wy7r). 


(4) 


A generic expression for the pion flux factor can be written as follows: 

f«/p( X L,t) = i-%2(l - gj . )^(0)-2a T (t) ^ ^ ^ F 2 (t, X L ), (5) 

2 A correction accounting for the small, but non-zero Q 2 values is applied, based on the Monte Carlo generator 
information, as explained in section 1331 
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where a_p(0) is the Pomeron intercept, a^{t) = a' n (t — ml) is the pion trajectory, g 2 pjTn /47 t is 
the pun coupling constant known from phenomenological analysis If25l of low energy data, and 
F(t, xl) is a form factor accounting for off mass-shell corrections and normalised to unity at 
the pion pole, F(ml, xl) = 1. There exists a variety of models for the exact form of the pion 
flux Il26l432ll which typically leads to a ~ 30% spread in the predicted cross section according 
to equation ©. Most of models use a non-Reggeized version of equation ©, i.e. ajp(0) = 1 
and a n (t) = 0. 


3 Experimental Procedure and Data Analysis 

3.1 HI detector 

A detailed description of the HI detector can be found elsewhere ll33l . Only those components 
relevant for the present analysis are described here. The origin of the right-handed HI coordi¬ 
nate system is the nominal ep interaction point. The direction of the proton beam defines the 
positive z-axis; the polar angle 6 is measured with respect to this axis. Transverse momenta 
are measured in the x-y plane. The pseudorapidity is defined by p = — In [tan(0/2)] and is 
measured in the laboratory frame. 

The central region of the detector is equipped with a tracking system. It included a set of 
two large coaxial cylindrical drift chambers (CJC), interleaved by a z chamber, and the central 
silicon tracker (CST) ll34l operated in a solenoidal magnetic field of 1.16 T. This provides 
a measurement of the transverse momentum of charged particles with resolution <j(pt)/pt — 
0.002 pt® 0.015 (pr measured in GeV), for particles emitted from the nominal interaction point 
with polar angle 20° < 9 < 160°. The interaction vertex is reconstructed from the tracks. The 
five central inner proportional chambers (CIP) ||35ll are located between the inner CJC and the 
CST. The CIP has an angular acceptance in the range 10° < 6 < 170°. The forward tracking 
detector is used to supplement track reconstruction in the region 7° < 9 < 30° and improves 
the hadronic final state reconstruction of forward going low momentum particles. 

The tracking system is surrounded by a finely segmented liquid argon (LAr) calorime¬ 
ter, which covers the polar angle range 4° < 6 < 154° with full azimuthal acceptance. The 
LAr calorimeter is used to measure the scattered electron and to reconstruct the energy of the 
hadronic final state. The backward region (153° <6 < 177.8°) is covered by a lead/scintillating- 
fibre calorimeter (SpaCal) |[36ll : its main purpose is the detection of scattered positrons. 

A set of “forward detectors” is sensitive to the energy flow close to the outgoing proton beam 
direction. It consists of the forward muon detector (FMD), the Plug calorimeter and the forward 
tagging system (FTS). The lead-scintillator Plug calorimeter enables energy measurements to 
be made in the pseudorapidity range 3.5 < p < 5.5. It is positioned around the beam-pipe 
at z = 4.9 m. The FMD is a system of six drift chambers which are grouped into two three- 
layer sections separated by a toroidal magnet. Although the nominal coverage of the FMD is 
1.9 < p < 3.7, particles with pseudorapidity up to 77 ~ 6.5 can be detected indirectly through 
their interactions with the beam transport system and detector support structures. The very 
forward region, 6.0 < p < 7.5, is covered by an FTS station which is used in this analysis. It 
consists of scintillator detectors surrounding the beam pipe at z = 28 m. The forward detectors 
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together with the LAr calorimeter are used here to suppress inelastic and proton dissociative 
background by requiring a large rapidity gap (LRG) void of activity between the leading neutron 
and the pions from the p° decay. 

Neutral particles produced at very small polar angles can be detected in the forward neutron 
calorimeter (FNC) B2lj37], which is situated at 106 m from the interaction point. It covers the 
pseudorapidity range rj > 7.9. The FNC is a lead-scintillator sandwich calorimeter. It consists 
of two longitudinal sections: the Preshower Calorimeter with a length corresponding to about 
60 radiation lengths, or 1.6 hadronic interaction lengths A, and the Main Calorimeter with a 
total length of 8.9A (see figure [2fr). The acceptance of the FNC is defined by the aperture of the 
HERA beam-line magnets and is limited to scattering angles of 9 £ 0.8 mrad with approximately 
30% azimuthal coverage, as illustrated in figure |2jx 

The absolute electromagnetic and hadronic energy scales of the FNC are known to 5% and 
2% precision, respectively 0. The energy resolution of the FNC calorimeter for electromag¬ 
netic showers is a(E)/E « 20%/y/E [GeV] © 2% and for hadronic showers a(E)/E pz 
63%/^E [GeV] © 3%, as determined in test beam measurements. The spatial resolution is 
cr{x,y) ph lOcm/i/E [GeV] © 0.6 cm for hadronic showers starting in the Main Calorimeter. 
A better spatial resolution of about 2 mm is achieved for electromagnetic showers and for those 
hadronic showers which start in the Preshower Calorimeter. 

The instantaneous luminosity is monitored based on the rate of the Bethe-Heitler process 
ep —y epp. The final state photon is detected in the photon detector located close to the 
beampipe at z = — 103 m. The precision of the integrated luminosity measurement is improved 
in a dedicated analysis of the elastic QED Compton process ll38ll in which both the scattered 
electron and the photon are detected in the SpaCal. 


3.2 Event selection 

The data sample of this analysis has been collected using a special low multiplicity trigger 
requiring two tracks with pr > 160 MeV and originating from the nominal event vertex, and 
at most one extra track with pt > 100 MeV. The tracks are found by the Fast Track Trigger 
(FTT) lf39l . based on hit information provided by the CJCs. The trigger also contains a veto 
condition against non-ep background provided by the CIP. The average trigger efficiency is 
about 75% for the analysis phase space. The trigger simulation has been verified and tuned to 
the data using an independently triggered data sample. 

For the analysis, exclusive events are selected, containing two oppositely charged pion 
candidates in the central tracker, a leading neutron in the FNC and nothing else above noise 
level in the dctcctoi[]. The photoproduction regime is ensured by the absence of a high en¬ 
ergy electromagnetic cluster consistent with a signal from a scattered beam positron in the 
calorimeters. This limits the photon virtuality to Q 2 < 2 GeV 2 , resulting in a mean value of 
(Q 2 ) = 0.04 GeV 2 . 

The p° candidate selection requires the reconstruction of the trajectories of two, and only 
two, oppositely charged particles in the central tracking detector. They must originate from a 

3 According to simulation, the forward going ir + from the proton vertex is emitted in the range q > 5.7 where 
it cannot be reliably measured or identified with the available apparatus. 
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common vertex lying within ±30 cm in z of the nominal ep interaction point, and must have 
transverse momenta above 0.2 GeV and polar angles within the interval 20° < 6 < 160°. The 
momentum of the p° meson is calculated as the vector sum of the two charged particle momenta. 
The two-pion invariant mass is required to be within the interval 0.3 < M nn < 1.5 GeV. Since 
no explicit hadron identification is used, events are discarded with Mkk < 1-04 GeV where 
M kk is the invariant mass of two particles under the kaon mass hypothesis. This cut suppresses 
a possible background from exclusive production of <f> mesons. 

Events containing a leading neutron are selected by requiring a hadronic cluster in the FNC 
with an energy above 120 GeV and a polar angle below 0.75 rnrad. The cut on polar angle, 
defined by the geometrical acceptance of the FNC, restricts the neutron transverse momenta to 
the range p T , n < %l ■ 0.69 GeV. 

To ensure exclusivity, additional cuts are applied on the calorimetric energy and on the 
response of the forward detectors. There should be no cluster with energy above 400 MeV, 
unless associated with p° decay products, in the SpaCal and FAr calorimeters. A Farge Rapidity 
Gap signature is required, by selecting events with no activity above noise levels in the forward 
detectors. This suppresses non-diffractive interactions to a negligible level and also significantly 
reduces diffractive background. 


Event selection (2006 — 2007, e + p) 

Analysis PS 

Measurement PS 

Trigger s 14 (low multiplicity) 

No e! in the detector 

Q 2 < 2 GeV 2 

Q 2 = 0 GeV 2 

2 tracks, net charge = 0, 

Pt> 0.2 GeV, 20° <6* <160°, 

20 < W 1P < 100 GeV 

20 < W 1P < 100 GeV 

from z vx < 30 cm 

Pt, p < 1.0 GeV 

-t' < 1.0 GeV 2 

0.3 < M nn < 1.5 GeV 

0.6 < M nn < 1.1 GeV 

2771^ < Mp < Mp + 5Tp 

ERG requirement 

~ 637,000 events 


E n > 120 GeV 

xl > 0.2 

0.35 < x L < 0.95 

6 n < 0.75 mrad 

6 n < 0.75 mrad 

PT,n < xl ■ 0.69 GeV 

~ 7000 events 

~ 6100 events 

~ 5770 events 

(OPE dominated range) OPE1 

Pt,u <0.2 GeV 

(~ 3600 events) 

OPE2 

PT,n<0.2GeV, 0.65<££<0.95 (~ 2200 events) 


Table 1: Event selection criteria and the definition of the kinematic phase space (PS) of the 
measurements. The measured cross sections are determined at Q 2 = 0 using the effective 
flux ©, based on the VDM. 

After these cuts the data sample contains about 7000 events. The event selection criteria 
together with the analysis and the measurement phase space definitions are summarised in ta¬ 
ble CD In order to better control migration effects and backgrounds most of the selection cuts 
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are kept softer than the final measurement phase space limits. In the end, the 7 p cross sections 
measured in the 6 n < 0.75 mrad range are based on ~ 5770 events. For the 771 cross section 
extraction additional cuts are applied in order to stay within a range where the validity of OPE 
can be safely expected. Two sub-samples are defined: OPE1 withp^n < 200 MeV, containing 
~3600 events and OPE2 withp T ,n < 200 MeV and xl > 0.65, containing ~2200 eventfl 

3.3 Monte Carlo simulations and corrections to the data 

Monte Carlo (MC) simulations are used to calculate acceptances and efficiencies for triggering, 
track reconstruction, event selection and background contributions and to account for migra¬ 
tions between measurement bins due to the finite detector resolution. 

Signal events from the DPP reaction (figure©) are modelled by the two-step MC generator 
POMPYT lf40l . in which the virtual pion is produced at the proton vertex according to one of 
the available pion flux parametrisations. This pion then scatters elastically on the photon from 
the electron beam, thus producing a vector meson (p° in our case). In this analysis the non- 
Reggeized pion flux factor is taken from the light-cone representation pT| with the form factor 
in equation © 

F 2 (t , x L ) = exp (^-Rln^z ’ ( 6 ) 

where R^n = 0.93 GeV~ 1 is the radius of the pion-proton Fock state 11271 . The same version 
of the pion flux factor has been used in previous HI publications on leading neutron measure¬ 
ments ©0 providing a good description of inclusive neutron spectra. For the numerical value 
of the pirn coupling constant, the most recent estimate If25l g 2 vw jAn r = 14.11 ± 0.20 is used. 

Since the exact shape of the p\ p dependence is not a priori known, two extreme versions 
are generated. In the first version a simple exponential shape is assumed, as expected for elastic 
p° photoproduction on the pion, with the slope 6 = 5 GeV -2 . For the second version a mass- 
dependent slope is taken, 4 < b(M n7r ) < 22 GeV -2 , typical for DPP processes as observed at 
lower energies IfTOlITT . 17]. The difference in the correction factors obtained using these two 
versions of MC simulations is part of the model dependent systematic uncertainty. 

The background events originating from diffractive p° production (figure ©l) are generated 
using the program DIFFVM Il42l . which is based on Regge theory and the Vector Dominance 
Model. All channels (elastic, single- and double-dissociation processes) are included, with 
the relative composition as measured in If43l . For the proton dissociative case the M Y mass 
spectrum is parametrised as da/dMp oc 1/M 2 - 16 , for My > 3.6 GeV 2 with quark and diquark 
fragmentation using the JETSET program PHI . For the low mass dissociation the production 
of excited nucleon states at the proton vertex is taken into account explicitly. Signal events, 
corresponding to the diagram shown in figure©, are excluded from the generated background 
sample. 

The DIFFVM program is also used to estimate possible contaminations from diffractive 
u;(782), 0(1020) and p'(1450 — 1700) production. 

As discussed in section [T] the pion exchange diagram dominates the cross section in the 
low t region where the contributions from the diagrams in figures © and |T]e almost cancel. 

4 The OPE2 sample corresponds to the low \t\ < 0.2 GeV 2 region, see figure[2j:. 
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To check a possible influence of these terms on the MC correction factors, neutron exchange 
events (b) were generated using POMPYT and events of class (c) using DIFFVM. As expected, 
these events have kinematic distributions and selection efficiencies similar to those from the 
pion exchange process and do not alter the MC correction factors beyond the quoted systematic 
uncertainties. 

In both the POMPYT and the DIFFVM generators a simple non-relativistic Breit-Wigner 
shape is used for the p meson mass. Therefore all MC events are reweighted to the relativistic 
Breit-Wigner shape with additional p T -dependent distortion as observed in p° photoproduction 
experiments. The distortion is caused by the interference between the resonant and non-resonant 
7 r + 7r _ production and is characterised by the phenomenological skewing parameter, urs, as 
suggested by Ross and Stodolsky l45l : 



with M p being the nominal resonance mass Pi6l and urs(pt,p) taken from published ZEUS data 
on elastic photoproduction of p° mesons 11471 . Additionally, the signal MC events (POMPYT) 
are reweighted in W 1P and in p\ to the observed shapes of the corresponding distributions. This 
reweighting is performed iteratively and has converged after two iterations. The uncertainty in 
the reweighting procedure is then taken into account in the systematic error analysis. 

Small, but non-zero values of Q 2 cause |£'| to differ from pj, by less than Q 2 . To account 
for this effect a multiplicative correction factor determined with the Monte Carlo generators is 
applied to the bins of the p\ distribution; the correction is obtained by taking the ratio between 
the |t'| and pj, distributions at the generator level. This correction varies from 1.1 at Pt, p = 0 
to 0.77 at p\ p = 1 GeV 2 . 

For all MC samples detector effects are simulated in detail with the GEANT program ll48l . 
The MC description of the detector response, including trigger efficiencies, is adjusted using 
comparisons with independent data. Beam-induced backgrounds are taken into account by 
overlaying the simulated events with randomly triggered real events. The simulated MC events 
are passed through the same reconstruction and analysis chain as is used for the data. 

The MC simulations are used to correct the distributions at the level of reconstructed par¬ 
ticles back to the hadron level on a bin-by-bin basis. The size of the correction factors is 12 
in average, corresponding to an efficiency of ~ 8%, and varies between ~ 10 and ~ 24 for 
different parts of the covered phase space. The main contributions to the inefficiency are: the 
azimuthal acceptance of the FNC (~ 30% on average), the p meson reconstruction efficiency 
which is zero if one of the tracks has low transverse momentum (~ 60%), the LRG selection 
efficiency (~ 60%) and the trigger efficiency (~ 75%). The bin purity, defined as the fraction 
of events reconstructed in a particular bin that originate from the same bin on hadron level, 
varies between 70% and 95% for one-dimensional distributions and between 45% and 65% for 
two-dimensional ones. As an example, figure |2fc illustrates the binning scheme used in the 
two-dimensional ( XL,PT,n ) distribution. 

3.4 Extraction of the p° signal 

The invariant mass distribution of the two tracks under the charged pion mass hypothesis is 
shown in figure [3Jl. The distribution is corrected for the mass dependent detector efficiency. 
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A fit is performed in the range M n7r > 0.4 GeV using the Ross-Stodolsky parametrisa- 
tion © for the p° meson mass shape and adding the contributions for the reflection from 
u —>• 7 T + 7 r - 7 r 0 and for the non-resonant background. Other sources of non-p° background, such 
as cu(782) —> 7 r + 7 r _ , 0(1020) —> K®Kg, 7 r + 7 r _ 7 r°, p' —> pnn, Air, titi, which may be misiden- 
tified as p° candidates, are estimated using MC simulations with the relative yield normalisa¬ 
tion fixed to previously measured and published values: a 7 P (o;)/< 7 7 p(p 0 ) = 0.10(±20%) l|49l , 
a 7 P (0)/a 7 P (p°) = 0.07(±20%) [j50]| and a 7 P (p')/cx 7 P (p 0 ) = 0.20(±50%) EJ. The resulting 
overall background contamination in the analysis region 0.6 < M <1.1 GeY is found to be 
(1.5 ±0.7)%. 

The fitted values of the resonance mass and width are 764±3(stat.) MeV and 155±5(stat.) 
MeV, respectively, in agreement with the nominal PDG values of M p and Y fl [|46l . The cross 
section is then calculated for the full mass range 2171^ < M nn < M p ± 5T P using the resonant 
part only, represented by the relativistic Breit-Wigner function BW p (M n7T ) with momentum 
dependent width T(M 7r7r ) Il52l : 


BW P (M. K7T ) 


M nn M p T(M 7r7r ) 

(M2-M2J2 + M 2 r(M _)2’ 


r(M 7r7r ) = v p 



(8) 


where q* is the momentum of the decay pions in the rest frame of a pair of pions with mass 
M n7r , and % is the value of q* for M n7T = M p . 

The Breit-Wigner shape is strongly distorted due to interference with the non-resonant tht 
production amplitude (dashed curve in figure©). The strength of the distortion is 7/7 -depcndcnt 
and within the ansatz © is characterised by the phenomenological skewing parameter, urs- 
For the full pr range of the present analysis, pip p < 1 GeV 2 , a fit results in the value n rs = 
4.22 ± 0.28. To study its pr dependence the fit is repeated in four p\ p bins. The values obtained 
are shown in figure©) in comparison with previously published ZEUS results ll47ll from elastic 
p° photoproduction, pp —)■ p°p. The dashed curve represents a fit to all these data by the 
empirical formula 

nRS = n 0 (p 2 T + M 2 ) _/3 (9) 

with no, M and /3 as free parameters. The fitted value of M 2 ~ 0.6 GeV 2 suggests that the 
relevant scale for photoproduction of vector mesons is indeed (py ± My). 

An important set of observables which characterise the helicity structure of the vector meson 
production are the angular distributions of the decay pions. Here we study the distribution of 
9 h which gives access to the p° spin-density matrix element r))). The angle 9 h is defined as the 
polar angle of the positively charged decay pion in the p° rest frame with respect to the meson 
direction in the centre-of-mass frame. According to the formalism presented in ll53ll the 
distribution 9h is given by: 


--— cx 1 - ± (3r{£ - 1) cos 2 9 h . (10) 

(J Cl COS u fa 

Figure © shows the acceptance corrected cos 9h distribution together with the fit by equa¬ 
tion (fTOl) yielding the value of Tqq = 0.108 ± 0.017. In figure© this result is compared to 
the values obtained in diffractive p° photo- and electro-production at HERA ll47l[540 551. The 
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steep Q 2 dependence is driven mainly by the QED gauge invariance motivated factor Q 2 /My, 
and can be fitted by a simple expression [|56l 


00 1 + £(M?/Q*)* 

with the parameters £ = 1.85 ± 0.10 and k = 0.67 ± 0.03, as illustrated by the dashed curve. 

In summary, all properties of the selected 7r + 7r - sample investigated here are consistent with 
p° photoproduction. 


3.5 Signal and background decomposition 

The event selection described in scction fOl does not completely suppress non-DPP background. 
According to the MC simulations, the remaining part is mostly due to proton dissociation with 
some admixture of double dissociative events. 

As in the case of inclusive leading neutron production Q signal and background events have 
different shapes of the leading neutron energy distribution, although in the present analysis the 
difference is less pronounced. The shape differences in the neutron energy spectrum predicted 
by MC for the DPP events (POMPYT) and for the proton dissociative background (DIFFVM) 
are still sufficient to disentangle these two contributions on a statistical basis. For this purpose a 
combination of the spectra obtained for reconstructed events of these two MC models fulfilling 
all selection criteria with free normalisation is fitted to the data. From this fit the background 
fraction is determined to be Fbg ~ 0 .34 ± 0.05. The uncertainty includes both the fit error 
and systematic uncertainties related to the background shape variation in terms of My and t 
dependencies and proton dissociation fraction in the overall diffractive cross section. Figure 0] 
illustrates this decomposition using the nominal DIFFVM parameters. 

Control plots for the data description by the Monte Carlo models using this signal to back¬ 
ground ratio are shown in figure [5} Since neither POMPYT nor DIFFVM are able to provide 
reliable absolute cross section predictions for such a final state, only a shape comparison is pos¬ 
sible. The irregular shape of the azimuthal angle distribution, ip n , is due to the FNC aperture 
limitations, as shown in figure [2b. 

In the fit described above the absolute normalisation for the DIFFVM prediction is left free. 
As a cross check, this normalisation has been fixed using an orthogonal, background dominated 
sample, obtained by requiring an ‘anti-FRG’ selection, i.e. p° +n events with additional activ¬ 
ity in the forward detectors. In this sample the background fraction is found to be 0.58 ± 0.07. 
Fixing the DIFFVM normalisation by a fit to the ‘anti-FRG’ sample results in a background 
contribution of F bg = 0.29 ± 0.05 in the main sample. Since the signal-to-background de¬ 
composition fit in this cross check gives a worse x 2 , the nominal value Fbg = 0.34 is used for 
the cross section determinations. The difference to the Fbg value determined in the nominal 
analysis, as described above, is well covered by systematic uncertainty of the FRG condition 
efficiency. 
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3.6 Cross section determination and systematic uncertainties 


The cross sections are measured for the kinematic ranges as defined in the rightmost column of 
table [U From the observed number of ep events, Ndata, the bin-integrated 7 p cross section in 
bin i is calculated as 



1 N l data ■ 

% £(A-e)i p 


( 11 ) 


where is the expected diffractive dissociation background in bin i taking into account the 
overall normalisation fraction Fbg ~ 0 .34, A ■ e is the correction for detector acceptance and 
efficiency, C is the integrated luminosity of the data, C p is the extrapolation factor for the num¬ 
ber of p° events from the M n7r measurement interval to the full p° mass range and <f > 7 = 0.1543 
is the value of the equivalent photon flux from equation © for the given (W rp , Q 2 ) ranged 
Since the statistics available does not allow for a reliable p° mass fit in every measurement bin, 
C l p is calculated using C p = 1.155, obtained from the fit of the full sample and bin-dependent 
skewing correction factor derived from the fitted dependence of nRs(pr , P ) ’ n equation ©. 


Several sources of experimental uncertainties are considered and their effects on the mea¬ 
sured cross section are quantified. The systematic uncertainties on the cross section measure¬ 
ments are determined using MC simulations, by propagating the corresponding uncertainty 
through the full analysis chain. The individual systematic uncertainties are grouped into four 
categories below. 


• Detector related sources. 

The trigger efficiency is verified and tuned with the precision of 3.4% using an indepen¬ 
dent monitoring sample. It is treated as correlated between different bins. 

The uncertainty due to the track finding and reconstruction efficiency in the central tracker 
is estimated to be 1% per track 11581 resulting in 2 % uncertainty in the cross section, taken 
to be correlated between bins. 

Several sources of uncertainties related to the measurement of the forward neutrons are 
considered. The uncertainty in the neutron detection efficiency which affects the mea¬ 
surement in a global way is 2% 0. The 2% uncertainty on the absolute hadronic energy 
scale of the FNC [Q leads to a systematic error of 1.1% for the xl -integrated cross sec¬ 
tion and varying between 2% and 19% in different xp bins. The acceptance of the FNC 
calorimeter is defined by the interaction point and the geometry of the HERA magnets 
and is determined using MC simulations. The uncertainty of the impact position of the 
particle on the FNC, due to beam inclination and the uncertainty on the FNC position, is 
estimated to be 5 mm Q. This results in an average uncertainty on the FNC acceptance 
determination of 4.5% reaching up to 10% for the pT, n distribution. 

The systematics due to the exclusivity condition in the main part of the HI detector is 
estimated to be 2.1%. It gets contributions from varying the LAr calorimeter noise cut 
between 400 MeV and 800 MeV (0.9%) and from the parameters of the algorithm con¬ 
necting clusters with tracks (1.9%). This error influences only the overall normalisation. 

5 Note, that the effective VDM flux (0 converts the ep cross section into a real 7 p cross section at Q 2 = 0, 
contrary to the EPA flux (57) converting it to the transverse j*p cross section, averaged over the measured Q 2 
range. The difference between the two approaches amounts to « 6 % integrated over the (Q 2 ■ y) range of the 
measurement. 
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The uncertainty from the LRG condition is determined in the same manner as in the 
HI inclusive diffraction analyses based on the large rapidity gap technique [f59l l60l. It 
is further verified by comparing the cross sections obtained using different components 
of the forward detector apparatus for the LRG selection: FMD alone vs FMD+FTS vs 
FMD+Plug vs FMD+FTS+Plug. The resulting uncertainty is conservatively estimated 
to be 9.0% affecting all bins in a correlated manner. 

• Backgrounds. 

Three different types of background are considered. 

Non-ep background is estimated from the shape of the ^-vertex distribution and from the 
analysis of non-colliding proton bunches to be (1.2±0.7)%. Background originating from 
random coincidences between p° photoproduction events and neutrons from p-gas inter¬ 
actions amounts to (1.0 ±0.2)%. This results in 2.2% background which was statistically 
subtracted in all distributions with an uncertainty of 0 . 8 %. 

Non-p° background, as discussed in section 13741 has an uncertainty of 0.7% and affects 
the overall normalisation only. 

Diffractive background to the DPP signal events (section l33T) is estimated with a precision 
of 7.6%. This is one of the largest individual uncertainties in the analysis. It is correlated 
between the bins. 

• MC model uncertainties. 

The uncertainty in the subtracted diffractive background due to the limited knowledge 
on 7 p diffraction is evaluated by varying the My and t dependencies in the DIFFVM 
simulation and the relative composition of diffractive channels within the limits allowed 
by previous HERA measurements. The resulting uncertainty is a part of the background 
subtraction systematics listed above. 

The systematic uncertainty of the MC correction factors for signal events is 4.1%, varying 
between 1% and 9% in different bins. It is evaluated from the difference between two 
versions of the POMPYT MC program with different p| dependencies of the p° cross 
section, as described in scction [±3l Here the uncertainty due to the POMPYT reweighting 
procedure is also accounted for. 

• Normalisation uncertainties. 

The uncertainty related to the p° mass fit, extrapolating from the measurement domain 
0.6 < M nn <1.1 GeV to the full mass range 2 m n < M n7r < M p ± 5T p , which implies a 
correction factor of C p = 1.155 on average in equation (fTTl) with an uncertainty of 1 . 6 % 
due to fit errors. 

The integrated luminosity of the data sample is known with 2.7% precision lf38l . 

Together with other normalisation errors listed above the resulting total normalisation 
uncertainty amounts to 4.4%. 

The systematic uncertainties shown in the figures and tables are calculated using the quadratic 
sum of all contributions, which may vary from point to point. They are larger than the statistical 
uncertainties in most of the measurement bins. 

The total systematic uncertainty for the integrated yp cross section is 14.6% including the 
global normalisation errors. 
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4 Results 


Total, single- and double-differential photoproduction cross sections for the reaction 7 p —> 
p°mr + are measured in the kinematic range defined in table Q] The photon-pion cross section, 
a 77T = cr(77r + —* p°7r + ) , is extracted from the differential cross section da 7P / d xl using the pion 
flux ll27ll integrated over the range pr, n < 0.2 GeV. The results are summarised in tables [2}(9] 
and are shown in figures l6lfT3l 


4.1 7 p cross sections 

The 7 p cross section integrated in the domain 0.35 < x j j < 0.95 and —t! < 1 GeV 2 and averaged 
over the energy range 20 < W 1P < 100 GeV is determined for two intervals of leading neutron 
transverse momentum: 


< 7(7 p —y p°nir + ) 

= (310 ± 6 sta t ± 45 sys ) nb 

for p T) n < xl ■ 0.69 GeV 

( 12 ) 

<j(pp -A p°n7T + ) 

= (130 ± 3 sta t ± 19 sys ) nb 

for pr,n < 0.2 GeV. 

(13) 


Single differential cross sections as a function of xl for these two regions are given in table |2] 
and are shown in figure [ 6 ] The data are compared in shape to the predictions based on different 
models for the pion flux. Some models, like FMS 11311 and NSSS Ii32l are disfavoured by 
the data and can be ruled out even on the basis of shape comparison alone. The other pion flux 
parametrisations: Bishari-0 ll26l . Holtmann |[27ll . KPP ||29l and MST lf30l are in good agreement 
with the data in both pj /n ranges. 

Additional constraints on the pion flux models could be provided by the dependence on t 
(or Pr n ) of the leading neutron. The double differential cross section d 2 a JP /dx L dpT . is mea¬ 
sured, and the results are presented in tableland figure [7] The bins are chosen such, that the 
data are not affected by the polar angle cut (see figure [2fc). Although neither the t-, nor the 
pj n -dependence of the pion flux models are exactly exponential they can be approximated by 
a simple exponent in many cases. Such an approximation has been used already in other analy¬ 
ses mm. The pf, ^distributions measured here for fixed xl are compatible with an exponential 
shape within the statistical and uncorrelated systematic errors. Therefore, the same approach 
is used here. The cross sections are fitted by a single exponential function er bnl - XL>p ' r - n in each 
xl bin. The quality of the fits is good, with P(x 2 ) — 0-35 A 0.60. The results are presented 
in table 0] and figure [ 8 j The measured 6 -slopes are compared to those obtained from several 
pion flux parametrisations. Despite of the large experimental uncertainties none of the mod¬ 
els is able to reproduce the datqj. A possible reason for this discrepancy could be the effect 
of energy-momentum conservation affecting the proton vertex in this exclusive reaction more 
strongly than in inclusive production of a leading neutron in which an apparent factorisation of 
the proton vertex has been observed. Another explanation ||T71[68 ;1 could be absorptive correc¬ 
tions which modify the t dependence of the amplitude, leading to an increase of the effective 
6 -slope at large xl as compared to the pure OPE model without absorption. 

6 Reweighting the signal MC using the measured b n (xL) slopes has only small effects on the cross section 
determination and is covered by the systematic uncertainties assigned to the pion flux models 


15 



The energy dependence of the reaction yp —* p°n7r + is presented in tableland in figure [9] 
The cross section drops with W 1P in contrast to the POMPYT MC expectation, where the energy 
dependence is driven by Pomeron exchange alone. A Regge motivated power law fit to the 
data, <J 1P (W 1P ) cc W^ p , yields 5 = —0.26 ± 0.06. s/ot ± 0.07 S3/S . The difference in the energy 
dependence in data and MC is also reflected in the pseudorapidity distribution of the p° meson, 
which is given in tableland shown in figure ITOl 

Finally, the cross section as a function of the four-momentum transfer squared of the p° 
meson, t', is given in table [7] and presented in figure HU It exhibits the very pronounced feature 
of a strongly changing slope between the low-t' and the high-t' regions. The fit is performed to 
the sum of two exponential functions: 

% = 0l e w ' + 02 e« (14) 

dr 

and yields the following slope parameters: 

b x = (25.72 ± ?>.22 unc ± 0.26 cor ) GeV -2 ; b 2 = (3.62 ± 0.30 unc ± 0.10 cor ) GeV -2 (15) 

where the first errors include statistical and uncorrelated systematic uncertainties and the second 
errors are due to correlated systematic uncertainties. In a geometric picture, the large value of b\ 
suggests that for a significant part of the data p° mesons are produced at large impact parameter 
values of order (r 2 ) = 2 by- {he) 2 ~ 2fm 2 « (1.6i7 p ) 2 . In other words, photons find pions 
in a cloud which extends far beyond the proton radius. The small value of b 2 corresponds to 
a target size of ~ 0.5 fm. In the DPP interpretation lfl7l[20l[2Tll the observed behaviour is a 
consequence of the interference between the amplitudes corresponding to the diagrams a, b and 
c in figure |U leading to a slope dependence on the invariant mass of the (nn + ) system produced 
at the proton vertex. Since the forward pion is not detected in this analysis the ( nn + ) invariant 
mass cannot be determined with sufficient precision, which prevents explicit measurement of 
the b(m) dependence. 

In order to investigate the presence of a possible factorisation between the proton and the 
photon vertices, the t' distribution is studied in bins of xl- The result of the fit by equation (fl4l) 
with x L dependent parameters a ; (.xyj, 5 ( (.xyJ is presented in tableland in figure HU in com¬ 
parison with the values given in equation (fT5l) for the full xl range. Also the evolution with 
x L of the ratio of two components, a\/a 2 , where a* = ^(1 — e~ bi ), is shown. Given the large 
experimental uncertainties no strong conclusion about factorisation of the two vertices can be 
drawn. 


4.2 77r cross section 

The pion flux models compatible with the data in shape of the xl distribution are used to ex¬ 
tract the photon-pion cross sections from da/d xl in the OPE approximation. The results are 
presented in tableland in figure H3 As a central value the Holtmann flux ll27ll is used, and the 
largest difference to the other three predictions lt26l|29l(30l provides an estimate of the model 
uncertainty which is ~ 19% on average. From the total yp cross section in equation (fT3l) and 
using the pion flux ([5]l6]) integrated in x L and p T , n , 1% = 0.056, the cross section of elastic 
photoproduction of p° on a pion target is determined at an average energy (W 77r ) ~ 24 GeV: 

a( 77 r + —»■ p° 7 r + ) = (2.33 ± 0.34(exp)7a 40 (model)) /ib, (16) 
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where the experimental uncertainty includes statistical, systematic and normalisation errors 
added in quadrature, while the model error is due to the uncertainty in the pion flux integral 
obtained for the different flux parametrisations compatible with our data. 

Theoretical studies of leading neutron production in ep collisions [[271 [29] suggest that in 
addition to the pion exchange process other processed may contribute at 10 — 20% level. To 
suppress these contributions it is recommended to perform cross section measurements in the 
‘OPE safe’ phase space region: low p Tn and high xl- In order to investigate a possible influence 
of non-OPE contributions the extraction of the photon-pion cross section is repeated for two 
additional regions, in which the validity of pure OPE is assumed. The cross sections for the 
full FNC acceptance range (9 n < 0.75 mrad, 0.35 < xl < 0.95) and for the OPE2 sub¬ 
sample (pT,n < 200 MeV, 0.65 < xp < 0.95) together with the value (fl6l) obtained for the 
OPE1 sample are presented in table[T0l The values of a(p7r + —* p°7r + ) extracted in these three 
different phase space regions agree well within the experimental errors. Thus no evidence for an 
extra contribution beyond the OPE is found in the full FNC acceptance range for the exclusive 
reaction studied here. 

Taking a value of cx(yp —>• p°p) = (9.5 ± 0.5) /ib at the corresponding energy (W) = 24 
GeV, which is an interpolation between fixed target and HERA measurements (see e.g. fig¬ 
ure 10 in (471), one obtains for the ratio r el = a^/cbf = 0.25 ± 0.06. A similar ratio, but for 
the total cross sections at (W) = 107 GeV, has been estimated by the ZEUS collaboration as 
r t ot = ^tot/ a tot = 0.32 ± 0.03 [El. Both ratios are significantly smaller than their respective ex¬ 
pectations, based on simple considerations. For r to t, a value of 2/3 is predicted by the additive 
quark model lf6TTl . while r e i = (^ 2E ) • (cr/^/cq^) 2 = 0.57 ± 0.03 can be deduced by combining 
the optical theorem, the eikonal approach (621 relating cross sections with elastic slope param¬ 
eters (631 and the data on pp, ir + p (641 and yp (471 elastic scattering. Such a suppression of 
the cross section is usually attributed to rescattering, or absorptive corrections (6514681 . which 
are essential for leading neutron production. For the exclusive reaction yp —>■ p°n7r + studied 
here this would imply an absorption factor of K abs = 0.44 ± 0.11. It is interesting to note, 
that this value is similar to the somewhat different, but conceptually related damping factor in 
diffractive dijet photoproduction, the rapidity gap survival probability, (S' 2 ) ~ 0.5, which has 
been determined by the HI collaboration [(6914711 . 


5 Summary 


The photoproduction cross section for exclusive p° production associated with a leading neutron 
is measured for the first time at HERA. The integrated yp cross section in the kinematic range 
20 < W 1P < 100 GeV, 0.35 < x L < 0.95 and 9 n < 0.75 mrad is determined with 2% statistical 
and 14.6% systematic precision. The elastic photon-pion cross section, cr(y-7r + —> p°7r + ), at 
(Wyn) = 24 GeV is extracted in the one-pion-exchange approximation. 

Single and double differential yp cross sections are measured. The differential cross section 
da/dt' shows a behaviour typical for exclusive double peripheral exchange processes. 

7 For inclusive leading neutron production, p. 02 trajectories should be considered, while for the exclusive reac¬ 
tion © the diagrams shown in figures [Th,c become important at larger t. 
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The differential cross sections for the leading neutron are sensitive to the pion flux models. 
While the shape of the xl distribution is well reproduced by most of the pion flux parametrisa- 
tions, the xl dependence of the pr slope of the leading neutron is not described by any of the 
existing models. This may indicate that the proton vertex factorisation hypothesis does not hold 
in exclusive photoproduction, e.g. due to large absorptive effects which are expected to play 
an essential role in soft peripheral processes. The estimated cross section ratio for the elastic 
photoproduction of p° mesons on the pion and on the proton, r ei = cr^/cTf 1 = 0.25 ± 0.06, 
suggests large absorption corrections, of the order of 60%, suppressing the rate of the studied 
reaction 7 p —> p 0 nn + . 
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(pT,i 

i < x L 

• 0.69 GeV) 


(PT,n < 0.2 GeV) 

x L 

da lp /dx L 

[H 

®stat 

[%] 

£unc 

u sys 

[%] 

£cor 
v sys 

[%] 

Stot 

m 

da lp /dx L 

[/ib] 

&stat 

[%] 

£unc 
v sys 

[%] 

£cor 
v sys 

[%] 

Stot 

[%] 

0.35 - 0.45 

0.213 

9.8 

10.6 

15.1 

20.9 

0.119 

11.2 

10.3 

15.2 

21.5 

0.45 - 0.55 

0.398 

7.0 

9.8 

15.4 

19.5 

0.164 

8.6 

7.5 

15.3 

19.1 

0.55-0.65 

0.530 

5.9 

7.2 

15.7 

18.2 

0.190 

7.6 

7.8 

15.4 

18.9 

0.65 - 0.75 

0.761 

4.1 

6.9 

12.8 

15.1 

0.274 

5.1 

9.5 

12.0 

16.2 

0.75 - 0.85 

0.806 

3.6 

5.0 

11.7 

13.2 

0.354 

4.1 

5.8 

10.7 

12.8 

0.85-0.95 

0.402 

5.4 

19.4 

12.8 

23.9 

0.204 

6.3 

15.0 

11.2 

19.7 


Table 2: Differential photoproduction cross sections d a lp / d xl for the exclusive process 7 p —> 
p°nir + in two regions of neutron transverse momentum and 20 < W 1P < 100 GeV. The statistical, 
uncorrelated and correlated systematic uncertainties, S sta t, and S c s ° r s respectively, are given 
together with the total uncertainty S to t, which does not include the global normalisation error of 

4.4%. 


xl range 

(xl) 

Pr,n ran § e 
[GeV 2 ] 

(PT,n) 

[GeV 2 ] 

d z a lp 

Ostat 

[%] 

ZlLTlC 

u sys 

[%] 

vcor 

°sys 

[%] 

Stot 

[%] 

dx L dp| n 

[/ib/GeV 2 ] 

0.35-0.50 

0.440 

0.00-0.01 

0.00499 

3.178 

13.9 

6.3 

14.8 

21.3 



0.01 - 0.03 

0.01998 

3.545 

12.1 

5.4 

12.7 

18.4 



0.03 - 0.06 

0.04495 

2.974 

13.7 

6.1 

12.7 

19.7 

0.50-0.65 

0.581 

0.00-0.01 

0.00492 

5.242 

10.5 

4.3 

14.0 

18.0 



0.01 - 0.03 

0.01969 

4.925 

8.6 

4.1 

12.9 

16.0 



0.03 - 0.06 

0.04429 

3.344 

11.7 

4.7 

13.9 

18.8 



0.06-0.12 

0.08719 

2.775 

11.2 

7.3 

13.7 

19.1 

0.65-0.80 

0.728 

0.00-0.01 

0.00489 

9.623 

6.3 

4.5 

11.4 

13.8 



0.01-0.03 

0.01957 

7.229 

5.5 

5.5 

12.0 

14.3 



0.03 - 0.06 

0.04403 

5.333 

7.3 

5.7 

12.2 

15.3 



0.06-0.12 

0.08617 

2.927 

8.4 

4.8 

13.7 

16.8 



0.12 - 0.20 

0.15324 

1.494 

14.7 

6.3 

17.9 

24.0 

0.80-0.95 

0.863 

0.00-0.01 

0.00484 

7.990 

7.6 

8.5 

11.2 

16.0 



0.01-0.03 

0.01935 

6.457 

5.7 

7.1 

10.9 

14.2 



0.03 - 0.06 

0.04354 

3.850 

7.9 

7.4 

12.3 

16.4 



0.06-0.12 

0.08425 

1.580 

11.3 

7.8 

15.7 

20.8 



0.12-0.30 

0.16558 

0.520 

14.1 

9.3 

18.7 

25.2 


Table 3: Double differential photoproduction cross sections d 2 a, w /dxfdpj n in the range 20 < 
W 1P < 100 GeV. The statistical, uncorrelated and correlated systematic uncertainties, S sta t, $sys 
and Sg° r s respectively, are given together with the total uncertainty S to t, which does not include 
the global normalisation error of 4.4%. 
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xl range 

(tl) 

b n [GeV 2 ] 

0.35 - 0.50 

0.440 

2.23 ±4.57 ± 2.10 

0.50 - 0.65 

0.581 

8.51 ± 1.74 ± 1.10 

0.65 - 0.80 

0.728 

13.17 ±0.90 ±0.65 

0.80 - 0.95 

0.863 

18.21 ±0.94 ± 1.05 


Table 4: The effective exponential slope, b n , obtained from the fit of double differential pho¬ 
toproduction cross sections (\xfApj n to a single exponential function in bins of xl- The 

first uncertainty represents the fit error from the statistical and uncorrelated systematic uncer¬ 
tainty and the second one is due to the correlated systematic uncertainty. 


W 1P [GeV] 

«T 7 

o’( 7 p —* p°mi + ) [nb] 

20-36 

0.06306 

343.7 ± 10.1 ±45.4 

36-52 

0.03578 

308.7 ± 12.3 ±43.5 

52-68 

0.02413 

294.2 ±15.8 ±45.2 

68-84 

0.01769 

260.0 ±23.1 ±44.9 

84 - 100 

0.01362 

214.5 ±50.2 ±45.0 


Table 5: Energy dependence of the exclusive photoproduction of a p° meson associated with a 
leading neutron, yp —> p°rm + . The first uncertainty is statistical and the second is systematic. 
The global normalisation uncertainty of 4.4% is not included. <h 7 is the integral of the photon 
flux 0]) in a given W 1P bin. 
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Vp 

da yp /drj 

[nb] 

Ostat 

[%] 

£unc 
vsys 

[%] 

Stcor 
v sys 

m 

Stot 

[%] 

[-5.0; -4.5) 

0.9 

68 . 

28. 

12 . 

75. 

[—4.5; -4.0) 

5.1 

27. 

18. 

11 . 

34. 

[—4.0; -3.5) 

8.8 

22 . 

11 . 

12 . 

27. 

[—3.5; —3.0) 

23.7 

14. 

6.1 

12 . 

20 . 

[—3.0; —2.5) 

44.0 

9.7 

4.0 

13. 

17. 

[—2.5; —2.0) 

45.2 

9.3 

3.1 

16. 

18. 

[-2.0;-1.5) 

47.5 

8.6 

3.5 

17. 

19. 

[-1.5;-1.0) 

48.2 

7.6 

2.9 

15. 

17. 

[—1.0; —0.5) 

45.9 

7.1 

5.9 

13. 

16. 

[—0.5; 0.0) 

38.9 

8.0 

3.2 

14. 

16. 

0.0;+0.5) 

46.2 

6.9 

5.7 

13. 

16. 

[+0.5; +1.0) 

52.1 

6.7 

7.1 

13. 

16. 

[+1.0; +1.5) 

63.8 

6.0 

5.4 

13. 

15. 

[+1.5; +2.0) 

86.2 

5.8 

4.4 

13. 

14. 

[+2.0;+2.5) 

39.8 

7.7 

3.1 

12 . 

15. 

[+2.5; +3.0) 

17.7 

11 . 

4.0 

12 . 

17. 

[+3.0; +3.5) 

7.8 

17. 

6.8 

12 . 

22 . 

[+3.5;+4.0) 

3.4 

26. 

11 . 

12 . 

30. 

[+4.0; +4.5) 

1.0 

55. 

21 . 

11 . 

60. 

[+4.5; +5.0) 

0.7 

64. 

33. 

11 . 

73. 


Table 6 : Differential photoproduction cross section da lv j dr/ for the exclusive process 7 p —> 
p°mr + as a function of the p° pseudorapidity in the kinematic range 0.35 < x L < 0.95, 9 n < 
0.75 mrad and 20 < W yp < 100 GeV. The statistical, uncorrelated and correlated systematic 
uncertainties, S stat , 8%™ and 8™ r s respectively, are given together with the total uncertainty 5 tot , 
which does not include the global normalisation error of 4.4%. 
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—t' range 
[GeV 2 ] 

(-?) 

[GeV 2 ] 

da rp /dt' 

[pb/GeV 2 ] 

Ostat 

[%] 

xunc 
v sys 

[%] 

vcor 
^sys 

[%] 

Stot 

m 

0.00-0.02 

0.0094 

2.771 

4.5 

2.5 

12.1 

13.2 

0.02 - 0.05 

0.0338 

1.821 

4.9 

1.7 

13.0 

14.0 

0.05-0.10 

0.0727 

0.996 

5.9 

1.3 

14.6 

15.8 

0.10-0.15 

0.1236 

0.600 

8.7 

1.0 

16.3 

18.5 

0.15-0.20 

0.1741 

0.402 

11.6 

2.9 

17.8 

21.4 

0.20-0.25 

0.2242 

0.343 

12.0 

3.7 

16.0 

20.3 

0.25 - 0.35 

0.2973 

0.279 

8.6 

5.1 

13.8 

17.0 

0.35 - 0.50 

0.4189 

0.178 

8.3 

6.4 

12.7 

16.4 

0.50 - 0.65 

0.5689 

0.104 

9.2 

7.8 

11.6 

16.8 

0.65-1.00 

0.7924 

0.037 

9.4 

18.7 

11.5 

23.9 


Table 7: Differential photoproduction cross section da rp /dt' for the exclusive process 

7 p —> p°mr + as a function of the p° four-momentum transfer squared, t', in the ki nematic 
range 0.35 <x L < 0.95, 6 n < 0.75 mrad and 20 < W lp < 100 GeV. The statistical, uncorrelated 
and correlated systematic uncertainties, S stat , and 8^°J S respectively, are given together with 
the total uncertainty 5 to t, which does not include the global normalisation error of 4.4%. 


xl range 

(xl) 

b x [GeV~ 2 ] 

b 2 [GeV" 2 ] 

cri /cr 2 

0.35-0.50 

0.440 

18.6 ±4.2 

2.54 ± 0.79 

1.501 ± 1.024 

0.50-0.65 

0.581 

26.0 ±5.5 

2.79 ± 0.43 

0.782 ±0.316 

0.65-0.80 

0.728 

28.1 ± 7.9 

4.24 ± 0.34 

0.244 ±0.091 

0.80-0.95 

0.863 

27.9 ±6.5 

4.42 ± 0.50 

0.394 ±0.142 

0.35-0.95 

0.686 

25.7 ±3.2 

3.62 ±0.32 

0.492 ±0.143 

Pr : n range [GeV 2 ] 

[GeV 2 ] 

bi [GeV" 2 ] 

b 2 [GeV" 2 ] 

cri/cr 2 

0.0-0.04 

0.015 

26.8 ±4.5 

4.07 ±0.34 

0.384 ±0.077 

0.04-0.30 

0.092 

26.6 ±4.4 

3.08 ± 0.46 

0.635 ±0.423 


Table 8: Exponential slopes, b\ and b 2 , and the ratio oi/o 2 , obtained from the components of 
fit (fl4l) to the differential cross section da rp /dt' in bins of xl and in bins of />f n . The errors 
represent the statistical and systematic uncertainties added in quadrature. 
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x L range 

r^(x L ) 

(14V) [GeV] 

Cr(77T + —> p°7T + ) [pb] 

0.35 

-0.45 

0.04407 

34.08 

2.71 ± 0.58l^ 

0.45 

-0.55 

0.07262 

31.11 

2.25 ± 0.43toil 

0.55 

-0.65 

0.10400 

27.83 

1.83 ± 0.35to!23 

0.65 

-0.75 

0.13154 

24.10 

2.09 ± 0.341°;| 

0.75 

-0.85 

0.13386 

19.68 

2.65 ± 0.34 

0.85 

-0.95 

0.07431 

13.91 

2.74 ± 0.54 


Table 9: Energy dependence of elastic p° photoproduction on the pion, 771 + —» p°7r + , extracted 
in the one-pion-exchange approximation using OPE1 sample. The first uncertainty represents 
the full experimental error and the second is the model error coming from the pion flux uncer¬ 
tainty (see text). T n (xrJ represents the value of the pion flux d5}(6l) integrated over the p T . n < 0.2 
GeV range, at a given x L . 


x L range 

Ptr [GeV] 

rv 

(W^) [GeV] 

<t( 77T + —^ p°7T + ) [pb] 

0.35-0.95 

0.35 - 0.95 

0.65 - 0.95 

xl ■ 0.69 

0.2 

0.2 

0.13815 

0.05604 

0.03397 

23.65 

23.65 

19.73 

2.25 ± 0.34tolo 

2.33 ± 0.34 toil 

2.45 ± 0.33tg;^ 


Table 10: Cross section of elastic p° photoproduction on the pion, 77r + —>■ p°7r + , extracted 
in the one-pion-exchange approximation using three different samples: full sample, OPE1 and 
OPE2. The first uncertainty represents the full experimental error and the second is the model 
error coming from the pion flux uncertainty (see text). T^ represents the value of the pion flux 
(|5]{6]) integrated over the corresponding (.17, p T , n ) range. 
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Figure 2: Sketch of the Forward Neutron Calorimeter (a), acceptance in the azimuthal plane (b) 
and the (x L) p\ n ) plane (c). The shaded area in figure (b) is the projected aperture limited by the 
proton beamline elements. The insert of figure (c) shows the acceptance in terms of the xl and t 
variables. The dark green area indicates the OPE2 region, p T n < 0.2 GeV and xl > 0.65. The 
curve in the main figure corresponds to the angular cut 6 n < 0.75 mrad, and the grid shows the 
binning scheme used for the double differential cross section measurement, d 2 a lp /dxLdp\ n . 
The dots are events from the preselection sample described in the leftmost column of table H] 
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Figure 3: The p° meson properties: (a) Mass distribution of the 7r + 7r - system for exclusive p° 
production with p^ < 1.0 GeV 2 associated with a leading neutron. The data points are cor¬ 
rected for the detector efficiency. The curves represent different components contributing to the 
measured distribution and the Breit-Wigner resonant part extracted from the fit to the data. The 
analysis region 0.6 < M n+7T ~ <1.1 GeV is indicated by vertical arrows, (b) Ross-Stodolsky 
skewing parameter, urs, as a function of [> 2 r of the 7r + 7r“ system. The values measured in this 
analysis are compared to previously obtained results for elastic photoproduction of p° mesons, 
IP P°P> by the ZEUS Collaboration, (c) Decay angular distribution of the n + in the helicity 
frame, (d) Spin-density matrix element, Tqq, as a function of Q 2 for diffractive p° photo- and 
electro-production. The curves on figures (b-d) represent the results of the fits discussed in the 
text. 


28 























p° with Forward Neutron 



EJGeV] 



Figure 4: Signal to background decomposition in the selected data sample. On the left panel 
the distribution of the measured neutron energy, E n , is shown together with the contributions 
from signal and background. On the right panel the x 2 dependence on the background fraction, 
Fb g , is shown. The shaded band represents the la uncertainty around the optimal fit value of 
the Fbg, taking into account statistical errors, FNC calibration systematics and the uncertainty 
in proton dissociation background shape. 
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Figure 5: Distributions of the reconstructed quantities pj, and 77 of the p° meson, xl, Pt an d 
p of the neutron and W 1P for data and Monte Carlo simulations normalised to the data. Data 
points are shown with statistical errors only. The shaded band indicates the uncertainty in the 
estimated background fraction. 
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X L X l 

Figure 6: Differential cross section da lp /dxL in the range 20 < W 1P < 100 GeV compared 
to the predictions based on different versions of the pion flux models. Top row: cross sections 
in the full FNC acceptance 6 n < 0.75 mrad. Bottom row: cross sections for the OPE1 range, 
PT,n < 200 MeV. Left-hand column: disfavoured versions of the pion fluxes, right-hand column: 
pion fluxes compatible with the data. The data points are shown with statistical (inner error bars) 
and total (outer error bars) uncertainties, excluding an overall normalisation error of 4.4%. All 
predictions are normalised to the data. 
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Figure 7: Double differential cross section d 2 a lp /dx L dp\ n of neutrons in the range 20 < 
W-y P < 100 GeV fitted with single exponential functions. The cross sections in different x bins 
j are scaled by the factor 3 J for better visibility. The binning scheme is shown in figure [2]:. The 
data points are shown with statistical (inner error bars) and total (outer error bars) uncertainties 
excluding an overall normalisation error of 4.4%. 
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Figure 8: The exponential slopes fitted through the dependence of the leading neutrons 
as a function of xl- The inner error bars represent statistical errors and the outer error bars 
are statistical and systematic errors added in quadrature. The data points are compared to the 
expectations of several parametrisations of the pion flux within the OPE model. 
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Figure 9: Cross section of the reaction 7 p —> p°mr + as a function of W 1P compared to the 
prediction from POMPYT MC program, which is normalised to the data. The dashed curve 
represents the Regge motivated fit a oc W s with 5 = —0.26±0.06 sta t ±0.07 sys . The data points 
are shown with statistical (inner error bars) and total uncertainties (outer error bars) excluding 
an overall normalisation error of 4.4%. 
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Figure 10: The differential cross section da rp /dr/ as a function of pseudorapidity of the p° 
meson compared to the prediction from the POMPYT MC program, which is normalised to the 
data. The inner error bars represent statistical errors and the outer error bars are total errors 
excluding an overall normalisation uncertainty of 4.4%. 
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Figure 11: Differential cross section da lp /dt' of p° mesons fitted with the sum of two expo¬ 
nential functions. The inner error bars represent statistical and uncorrelated systematic uncer¬ 
tainties added in quadrature and the outer error bars are total uncertainties, excluding an overall 
normalisation error of 4.4%. 
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Figure 12: The two exponential slopes, b\ and 62 , obtained by fitting the t' dependence of the p° 
mesons (left) and the relative contribution of the two exponents to the overall cross section of 
the reaction pp —> p 0 mr + (right) as a function of xl (top) and pj n (bottom). The error bars re¬ 
present statistical and systematic uncertainties added in quadrature. Horizontal lines with error 
bands show the corresponding average values for the full ranges of 0.35 < xl < 0.95 and 
0 < p 2 T)U < 0.3 GeV 2 . 
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Figure 13: Elastic cross section, a® 1 ^ = < 7 ( 77 ? + —y p°n + ), extracted in the one-pion-exchange 
approximation as a function of the photon-pion energy, W , y7r . The inner error bars represent the 
total experimental uncertainty and the outer error bars are experimental and model uncertainties 
added in quadrature, where the model error is due to pion flux uncertainties. The dark shaded 
band represents the average value for the full W- jlT range as given in equation (fl6l) . 


38 









